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Using a perturbative approach, the effects of the energy gap induced by the Aharonov-Bohm (AB) 
flux on the transport properties of defective metallic single-walled carbon nanotubes (MSWCNTs) 
are investigated. The electronic waves scattered back and forth by a pair of impurities give rise to 
Fabry- Perot oscillations which constitutes a coherent backscattering interference pattern (CBSIP). 
' It is shown that, the CBSIP is aperiodically modulated by applying a magnetic field parallel to 

the nanotube axis. In fact, the AB-fiux brings this CBSIP under control by an additional phase 
■ shift. As a consequence, the extrema as well as zeros of the CBSIP are located at the irrational 

~r fractions of the quantity $p = &/&o, where <E> is the flux piercing the nanotube cross section and 

$o = h/e is the magnetic quantum flux. Indeed, the spacing between two adjacent extrema in 
the magneto-differential conductance (MDC) profile is decreased with increasing the magnetic field. 
The faster and higher and slower and shorter variations is then obtained by metallic zigzag and 
armchair nanotubes, respectively. Such results propose that defective metallic nanotubes could be 
used as magneto-conductance switching devices based on the AB effect. 
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PACS numbers: 72.10.-d, 73.63.Fg, 72.15.Rn, 85.35.Ds 



I. INTRODUCTION 



Due to their quasi-one-dimensional structure and intriguing electronic properties, carbon nanotubes have been at- 
tracted an increasing amount of attentions pj. Carbon nanotubes are tubular nano-objects which can be thought of 
as graphenes wrapped onto a seamless cylinder. Depending sensitively on the wrapping vector, a nanotube may be 
either a one-dimensional (ID) metal with a finite density of states at the Fermi energy or a semiconductor with a 
gap. In special, for the sake of the ID nature of their electronic conduction bands near the Fermi energy, metallic 
single-walled nanotubes constitute a nearly perfect realization of ID quantum wires 0, d, 0, 

The investigation of quantum transport in carbon nanotubes is expected to have unprecedented potential applica- 
tions for developing nanoelectronic devices. They can be applied as conducting quantum wires 0,0], single-electron 
tunneling tranistors0, [1], field-effect transistors Q , and spin-electronic devices [Toj. Theoretical calculations based 
on the Landauer-Biitticker formalism [Til. [T3 | predict the conductance quantization for a perfect metallic nanotube 
for the case of ideal contacts. The maximum value of the conductance near the Fermi energy reaches 2Gq, where 
Go = 2e 2 Jh is the conductance quantum[l3j|. However, in contrast to the pristine nanotube, several theoretical 
works [H, [U, [13, E 01, HI and experimental evidences [HI [H, HI HI have shown that in the presence of dis- 
' orders coming from various sources like chemical impurities, topological defects, Stone and Wales [251]. and vacancies 
this quantized conductance of the nanotube does not follow the aforesaid results. Practically, these imperfections are 
unavoidable in manipulation nanotubes into devices and induce departure from ballistic transport, and yet preserve 
• i-H , quantum interference effects, which can be profoundly affected by magnetic fields. 

Owing to the decoherence, the quantum corrections to the classical conductance of a device are usually negligible 
in macroscopic systems at the room temperature. In mesoscopic systems at low temperatures, however, the quantum 
mechanical coherency becomes more important because the phase coherence length 1$ increases with decreasing tem- 
perature. When the coherence length 1^ exceeds the elastic mean free path l m , scattering on different impurities can 
interfere. Several quantum interference (QI) modifications are (1) the WL correction, which originates from pairs of 
time-reversed paths in a diffusive sample interfere constructively in the zero magnetic field. This interference enhances 
(reduces) the probability of electronic backscattering, decreasing (increasing) the conductance of the sample [26|: (2) 
the AB and Altshuler- Aronov- Spivak (AAS) oscillations. The AAS effect is actually the same WL correction em- 
bracing a magnetic field[27l.[28l[29j. As the magnetic field is increased, the AB phase eliminates the WL constructive 
interference, leading to a magneto-conductance; (3) universal conductance fluctuations (UCF), which means that the 



'Electronic address: mh-bagheri@cc.sbu.ac.ir. 



2 



conductance fluctuations are independent of the conductor details. 

Furthermore, one of the unique properties of carbon nanotubes is that their metallicity can be controlled by 
an external magnetic field applied parallel to the nanotube axis. This magnetic field gives rise to a periodic en- 
ergy gap at the charge neutrality point (CNP), where the bonding and antibonding bands are crossed. When 
the cross section of a nanotube is pierced by the magnetic field, the electronic wavefunctions acquire an addi- 
tional phase 27r<I>/<I>o. Thus, metallic nanotubes can be made semiconducting and vice versa. Over the past few 
years, remarkable efforts have been undertaken to evidence the effects of a magnetic field on the band structure of 

nanotubes HE m M, HE H HEM i3, H Hi, SO, SH, SI, IH • 

Following our previous paper (44[, in which a perturbative approach is well developed to include effects of the band 
structure and impurity on transport characteristics of metallic nanotubes, the current work concentrates on elucidat- 
ing influences of the AB-flux [45j on the differential conductance (DC). The motivation of this attempt is to study 
how the magnetic field dependence of the band-structure of the nanotube influences the DC. This may provide us 
the possibility of fabricating magneto-conductance switching devices based on the AB effect in defective metallic nan- 
otubes. It is shown that, for a couple of impurities the nanotube behaves like a Fabry-Perot electron resonator [HE HH 
and the CBSIP resulting from the Fabry-Perot oscillations is aperiodically modulated in the presence of the AB-flux. 
Aperiodicity means that no specific magnetic flux periodicity is found in the MDC profile. Further, extrema as well 
as zeros of the MDC are positioned at irrational fractions of the magnetic flux with a spacing which is decreased by 
increasing the magnetic field. 

The paper is established as follow. In Sec. [TTJ the model of Ref.[44j is developed to include the AB-flux. In Sec. 
IIII1 we discuss the CBSIP in the presence of the AB-flux both for a single and for a couple of impurities. 



II. THEORETICAL MODEL 



We address a defective MSWCNT in the presence of an axial electric and magnetic field. The full descriptions 
of the model in the absence of the magnetic field can be found in Ref. [44|. Here, we just add the AB-flux in its 
band-structure, so the Hamiltonian of the whole system is given by 

H(® P ) = Ht u be($p)+H sd + H imp . (1) 
In the above equation the first term, describing the kinetic energy of electrons for a perfect nanotube, is given bv[44ll46j 
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In the presence of an uniform magnetic field B parallel to the nanotube axis, the wrapping modes are modified 
according to q/r t — > q/r t + <& p /r t [31], so the magnetic field dependent band-structure £^ + q, f (fc) is|4| 



± < 1 + 4 cos 



7o 



—^-clcc I — [<Z + 'PpJ sinw + fccosw 



3 (I 



+4 cos z 



^3 fl, _ , . 

—^~ a cc I — [q + Qpl smu) + kcosuj 



(3) 



where operators C^(k) and C^(k) create and destroy electrons in the orbital with energy £^ + $ (k), respectively. 
The + and — signs correspond to the conduction and valence band, respectively. Good quantum numbers of electron 
states are (q,k) where q = l,...,N t /2 and k € (-tt/T,tt/T). The quantities N t , M = N t /2, T, r u a cc ~ 1.44A, 
and 70 — 3.0 eV are the number of carbon atoms in the nanotube unit cell, the number of graphene unit cells in a 
given nanotube unit cell, the length of the translation vector, the nanotube radius, the C-C bond length, and the 
nearest- neighbor overlap integral energy, respectively. Also, lj — -k /Q — where 9 is the chiral angel of the nanotube 
whose value for the armchair and zigzag nanotube is 7r/6 and 0, respectively. It is assumed that the on-site energy is 
zero and the Fermi energy remains unchanged at the CNP. In the zero magnetic field, all metallic linear bands cross 
the undoped Fermi level either degenerated at kF — (metallic zigzag) or separated at kp = ±2ir/3T (armchair) 
in the first Brillouin zone (FBZ). $ p equals &/&o, with $ = nr^B. When <E> p becomes an integer, the AB-flux is 
canceled by q. It means that the gap induced by the magnetic field oscillates periodically and can be obtained by the 
expression A s ($ p ) = 2min{| £^ +(S , p (k) |}(see Fig. QJi). Lu[42|] has shown that, for metallic nanotubes the energy gap 
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induced by an axial magnetic field is expressed by 
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where Ao = 7oa C c/?"t defines a characteristic energy associated with the nanotube. Note that the expression 
£f + q, (k) — ±70 sin[7r(g + $ p )/n] gives van Hove singularities (VHSs) positions. Further, for later calculations, we 
have exploited the corresponding Bloch's states of an isolated nanotube previously derived in Ref . . 

For considering the magneto-transport properties of the nanotube near the Fermi level, we adopt the light-cone 
approximation of the dispersion relation of Eq. ([3]) which provides us a simple formula of the s-th ID subband around 
k F . Thus, Eq. © reduces to[H[ 
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where s = 1, . . . , N t /2. For the lowest lying subband, with s = 1 around k — ±kp, the energy band gap in the absence 
of the magnetic field is zero. Using eq. {5} one obtains A g ($ p ) = 3Ao$ p . As the field strength increases the line 
through the Fermi energy at zero magnetic field is shifted away further from the CNPs thus given rise to an increasing 
energy-gap. It is also worth mentioning that, the quantity /i or b = ev F r t /2, with v F — 3~/oa cc /2h ~ 10 6 m/s, is the 
magnetic moment of an electron traveling in a loop of radius r t with velocity v F . Changes in the energy of electron 
states can be described by the interaction of this orbital magnetic moment with an axial magnetic field. A magnetic 
field parallel to the nanotube axis is predicted to shift the energy of these states by AE = —fx or b ■ B = ±ev F rtB jl = 
±3A $ p /2 (see Fig. IHb)). 

Furthermore, the second and third terms in equation (1) are, respectively, the Hamiltonian of non-interacting 
electrons under the external sourcedrain voltage V s d and the Hamiltonian of the interaction of electrons with impurities 
[13, El, El| like those presented in Q. Eventually, upon substituting q in Eq. (20) of Ref. [3 by q + <& p , we obtain 
the dimensionless form of the MDC at the zero temperature as follow 
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where G™ al 
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Gtm P + G imp - Also, w* +$p (k) = (l/K)d£f + $ (k)/dk is the electron velocity, and J^ q af3 {k, k') is a ma- 
trix for the impurity potential located at a position, namely, xe |44j. We have also assumed that the magnetic field does 
not affect the Fermi energy, i.e. £r(B) = £f(0). More importantly, the expression [sign[v^ +<i , {k)]sign[v^, + ^, (k')] — 1] 
controls the scattering event from the initial state to the final state via the sign of the electron velocity. It requires that 
only backward scattering events are possible in one-dimensional systems like nanotubes. The coherent backscattering 
(CBS) of the electron is an effect that describes the appearance of a backscattered peak when the electron traveling 
in a time-reversed path self-interferes constructively in the backscattered direction. It means that the electronic wave 
is weakly localized[5(| HJ, [52j ■ 

By obtaining the solutions of the energy-momentum conservation equation, i.e. £q + $ (k) = £gi + $ p (k + g) where g 
is the transferred momentum, we now evaluate Eq. ([6]) at some special k-points in the FBZ. Using Eq. |3]) for the 
(n, n) armchair nanotubes one obtains 
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and for the (n, 0) zigzag nanotubes the equivalent expression is given by 
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For the intrasubband scattering, i.e. \q + <& p , k) —> \q' + <5> p , fc') = \q + ^p, fc + g), Eq. ([7]) has four scattering roots as 
follow 
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while for the metallic zigzag nanotubes Eq. © provides only two roots and — 2k. 

k! = k : The root g ± = means that q and k are conserved, and no scattering event is occurred. Thus, the MDC 
becomes zero. 

k! = —k : The root g = — 2k describes the CBS of the electron within the same subband to another Fermi point. In the 
CBS effect, the electron is elastically scattered back to a momentum directly opposite to its original momentum state 
in the momentum space. Let later on replace Gi mp by Gqbs- For a couple of impurities located at — T/ x +Rj 1 +di 
and x n = f h + R J2 + d 2 , Eq. © yields [ill 



Gq 



, N t /2 FBZ 



8=1 fee I 



1 2 








2 


2 j 



x cos < 2k 



(l 2 -h)T+(R n -R n 



T 
T 



(10) 



where M. is the total number of nanotube unit cells[44|. Because (k) — —£ q+<s , (k); if £f(0) = then G,tg S = 

Gcbs- For the case of a single impurity the CBSIP is killed. Because two carbon atoms A and B inside a graphite 
unit cell belong to two different sublattices, the impurity can occupy one of the lattice site. For simplicity, we have 
here assumed that two impurities are substituted on B— sites with the same circumferential angle along the nanotube 
axis [44]]. These arrangements of impurities break all mirror symmetry planes containing the nanotube axis[l(|. By 
turning the sum over k into an integral and with exploiting Eq. §5§ for the lowest lying subband, Eq. (JTUJ) leads to 
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From Eq. (fTTj) one can draw several conclusions: (1) for 



the armchair nanotubes we have X = 1, Y = arm, and kp = 27r/3T a ™, while for the metallic zigzag ones X = 2, 
Y = zig, and kp = 0; (2) the cosine term is responsible for the CBSIP. Averaging over different impurity configurations 
melts away this interference term; (3) no switching effect from positive to negative MDC is occurred by changing 
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the orientation of the magnetic field with respect to the nanotube axis. This means that, the reciprocity relation 
GcBsi&p) — Gcbs{—$ p ) is fulfilled; (4) the amplitude of this CBSIP depends on both the source-drain voltage and 
the AB-flux; (5) in the limit <E> P — > 0, one recovers the solution of the free-magnetic field case derived in Ref.[4J]; (6) 
conduction through this gapped nanotube is dependent sensitively on the exact position of the V s d with respect to 
the lowest level subband edges. Strictly speaking, there is a threshold voltage determined by eV s d > 3Ao$ p /2 and 
eV s d < — 3Ao ( l ) p/2, below and above which, respectively, the transport is forbidden. This issue is in agreement with 
the density of state due to the one dimensional subbands expected for semiconductor nanotubes. In other words, the 
MDC is singular at the position of the lowest subband bottom indicating its van-Hove singularity; (7) a closer look 
at the argument of the second cosine term reveals that the interference term leads to aperiodic oscillations in the 
MDC profile. This is because this argument is a nonlinear mapping of the AB-flux as well as the source-drain voltage. 
In fact, the DC is aperiodically modulated through the AB-flux. At zero temperature, it would be plausible if we 
suppose that the system size plays the role of the phase-coherence length. In the presence of the AB-flux the electrons 
acquire additional phases, and we can control the interference pattern made from the conjugated time-reversed paths. 
More importantly is the negative MDC. Actually, it originates from not only the QI effects but also the pseudospin 
conservation rule. The negative MDC feature may be exploited for designing magneto-conductance switches based 
on the AB effect. 

k' = ±(2/^/3a cc ) arccos [cos(3(g + $ p )a cc /2r t ) + cos(v / 3fca C c/2)] : These two last roots are actually the intersubband 
backscattering around the the same Fermi point, and we currently discard them[56|. 



III. DISCUSSIONS 



Using the two-terminal Landauer-Butticker approach for a two-band model, the whole resistance of the nanotube 
is approximately given byJUJ 
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tube 
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In the above equation, the first term is the resistance of a perfect ballistic nanotube with perfect contacts. It originates 
from the redistribution of electrons between reservoirs and the nanotube. The second terms is the quantum correction 
coming from the CBS effect. Two last terms, discarded here, are for imperfect contacts between the nanotube and 
reservoirs. To investigate the behavior of the MDC as a function of the AB-flux, we have numerically performed 
Eq. (jlip for both armchair and zigzag nanotubes. Results arc the same for both repulsive and attractive impurity 
potentials. Let us suppose g=10 4 7o representing a typical impurity and M. — 1000000. In Eq. (fTTj) . the product of two 
cosine terms is actually a resultant wave coming from the superposition of two standing waves with the same amplitude 
but different wavenumbers k\ — k F + {l/2)^{eV sd /hv F ) 2 - (<& p /r t ) 2 and k 2 = k F - (l/2)^(eV sd /hv F ) 2 - (<5> p /r t ) 2 . 
These two initial standing waves describing two degenerate resonant states induced by impurities in the FBZ are 
given by 
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Because two functions fi and f 2 are not periodic in the $ p -space, so is their superposition, i.e. f\ + fi. Thus, an 
aperiodic variation in the MDC is expected. The phase difference for an electron propagating over the length l m is 
given by Sip(^ p ) = 2Akl m where Afc = ki — k 2 = \/{eV s d/fiv F ) 2 — ($ p /r t ) 2 . Constructive interference occurs when 
the extrema of two waves add together and the phase difference becomes an integer multiple of 7r, i.e. 5tp(& p ) = <jtt, 
with <7 6 Z. On the other hand, destructive interference occurs when two waves have a phase difference of a half- 
integer multiple of 7r, i.e. Sip(^ p ) = (a + l/2)ir. An analytic expression in the <E> p -space can be derived easily as 
follow 
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The above equation gives actually the spacing between the MDC extrema (constructive) or zeros (destructive)in the 
$ p -space. Due to the nonlinear mapping between Sip(^ p ) and $ p , the MDC vs. $ p behaves aperiodically. The most 
important feature is that, extrema and zeros are located at irrational fractions in the $ p -space. In other words, 
$ =irrationalx(&Q. The nonlinear dependence of the extrema positions as a function of a is depicted in Figs, dja) 
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and 2(b) for I2 — h = 50 and 250, respectively. In both panels the nonlinear behavior of oscillations can be seen by 
comparing the spacing between two horizontally adjacent lines. It should be pointed out that, variations in the MDC 
are aperiodic in the eV s d-sp&ce as well. 

Also, Ak can be expressed in terms of the series Afc = A (l + x/2 — x 2 /8 + x 3 /16 — •■•), where A = \eV s d/hvp\ 
and xi^p) = —($p/'''t) 2 /(eV s d/tiVF) 2 - Thus, the phase accumulated by the electron can be expressed by 
Sip = Sip(0) + 6<p($p), where 6<p(0) = 2Xl m is the phase difference in the absence of the AB-flux and 
S(p($> p ) = \l m (x — X 2 /4 + x 3 /8 — • • • ) is the magnetic field dependent phase difference. As a check, we see 
that, for x — 0, i.e. in the absence of the AB-flux, 5ip($ p ) = 0. For the (n, n) armchair and (n, 0) zigzag nanotubes 
we find r^ rm = 3na cc /27r; T arm — \/3a cc , and r\ 19 — \* r 3na cc /2Tr; T zlg = 3a cc , respectively, so in case we approximate 
6ip($ p ) by Xl m x, it is straight to show that 5(p($ p ) arm /5ip($p) zi9 = s/3. 

For a single impurity where 1% = l\, the nanotube is less defective. In this case the quantum interference due to 
the CBS is killed, and the AB modulation becomes dominant. Evaluating Eq. (TTTj) leads to a U-like behavior for the 
armchair and zigzag nanotubes. As depicted in Fig. [31 these curves are symmetric and centered at B = (or <5> p = 0). 
Each curve has a plateau which decreases with increasing the magnetic field and its magnitude depends strongly 
on the location of the source-drain voltage. A deep look at upon panels 3(a) and 3(b) shows that the magnitude of 
the zigzag plateau is approximately twice the armchair one for a fixed value of the source-drain voltage. Recently, 
Lassagne et al [1(3] have also observed such U-like curves, of course through a Schottky barrier for different gate 
voltages at non-zero temperatures, for a clean multi-walled nanotube threaded by the AB-flux. We should emphasize 
that, although our result share some similarities in the AB-pattern with that of Ref.[40j], but their underlying physical 
transport phenomena could be different. It is expected that, at a non-zero temperature and gate voltage such U-like 
behavior would drastically changed in our model. 

Moreover, for a couple of impurities, with l\ ^ I2, the MDC as a function of <& p for two different interdistances 
between impurities is calculated. In Fig. [4l panels show aperiodic fluctuations which alter between positive and 
negative values. The amplitude of oscillations is increased with increasing the magnetic field, while the spacing 
between two adjacent extrema is decreased. These fluctuations represent a hallmark of defective quantum transport 
resulting from the competition between the CBS effect and the AB-flux. Such fluctuations may be attributed to the 
Fabry-Perot oscillations [HJ modulated by the AB-flux. Positions of some extrema in the MDC are labeled by arrows 
in Fig. IHa). The most striking and immediately visible difference between armchair oscillations and zigzag ones, 
say by comparing panels 4(b) and 4(d), is that for the same value of the source-drain voltage the fluctuations of the 
zigzag nanotubes are faster and higher than that of the armchair one. The envelope functions of the extrema have a 
U-form portrait as well. 

In summary, this semi-classical study shows the subtle interplay between the quantum interference phenomena 
originating from Fabry-Perot oscillations 53] and the magnetic field-dependent of the band structure in defective 
metallic nanotubes. We have shown that, how such oscillations can be modulated using the AB-flux. Nonlinear 
mapping between the MDC and the magnetic filed leads to aperiodic fluctuations. Such results may be applied for 
manipulating defective metallic nanotubes into quantum interference devices, say, for the construction of nanotube 
magneto-conductance devices based on the AB effect 55J. Moreover, it is worth mentioning that the model is flexible 
to incorporate inelastic events like the electron-electron and electron-phonon scattering events. In the presence of 
such decoherent effects we expect a drastic change in the interference pattern of the differential conductance. 
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FIG. 1: (Color online) The energy dispersion relation for subbands of the armchair and metallic zigzag nanotubes in the presence 
of an axial magnetic field, (a) A multiband model, coming from evaluating Eq. <(3j in the presence of a 1000 Tesla magnetic 
field pointing along its axis, for the (6, 6) armchair (circle) and (6, 0) zigzag (square) nanotubes. The nanotube now has a finite 
subband-gap A 9 ($ p ) expressed by Eq. Q, and that all degenerate levels have been split. Antibonding bands (green- E q (k) > 0) 
are symmetric to the bonding bands (purple-_B q (fc) < 0). (b) A two-band model, which comes from evaluating Eq. ((5} (E s (k) 
with s = 1) in the presence of a 10 milli- Tesla magnetic field pointing along its axis, includes the (6, 6) armchair (solid line) and 
(6,0) zigzag (dotted) nanotubes. The subband-gap is now expressed by A 9 ($ p ) = 3Ao<3? p . The electron scattering processes 
change electrons from the right moving to the left moving leading to electrical resistance. Generally, both intrasubband and 
intersubband scattering events are likely. Energies are scaled in Rydberg and lengths in Bohr radius. 
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FIG. 2: (Color online) The positions of extrema in the <E> p -space are calculated with using the upper part of Eq. (|14|l . (a) 
The allowed cr's for the (6,6) (diamond-blue) and (6,0) (circle-red) nanotubes, with fa — h ~ 50 and eV a d = 0.11, are 18 and 
33, respectively, (b) The allowed cr's for the (6, 6) (diamond-blue) and (6, 0) (circle-red) nanotubes, with fa — h = 250 and 
eV s d = 0.11, are 93 and 160, respectively. The spacing between two horizontally adjacent lines is decreased with increasing the 
magnetic field which obviously shows that oscillations are aperiodic. 
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FIG. 3: (Color online) Calculated MDC as a function of the magnetic field B for a single impurity. Results come from evaluating 
Eq. for the armchair and metallic zigzag nanotubes. (a) Traces are plotted for the (6,6) armchair, with eV s d = 0.11 and 
B £ [-22,22] milli-Tesla (solid-blue), and eV sd = 0.15 with B G [-30,30] milli-Tesla (dotted-red). (b) Traces are plotted for 
the (6,0) zigzag, with eV sd = 0.11 and B G [-38,38] milli-Tesla (solid-blue), and eV ad = 0.15 with B G [-52,52] milli-Tesla 
(dotted-red). They exhibit a U-like behavior. The plateau of the zigzag nanotube is approximately twice the plateau of the 
armchair one. 
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FIG. 4: (Color online) The calculated MDC coming from evaluating Eq. for the nanotube (n, m) for a pair of impurities. 
The CBSIP shows aperiodic oscillations. The extrema as well as zeros are located at the irrational fractions of $ p . Positions of 
some extrema are indicated by arrows, (a) (n, m) = (6, 6), eV s d = 0.11, B £ [—22, 22] (mT), and h — h = 50; (b) (n, m) = (6, 6), 
eV sd = 0.11, B G [-22,22] (mT), and h - h = 250; (c) (n,m) = (6,0), eV sd = 0.11, B G [-38,38] (mT), and h~h = 50; (d) 
(n,m) = (6,0), eVsd = 0.11, B G [—38,38] (mT), and fa — h = 250. A comparison between, say panels (b) and (d), exhibits 
that the faster/higher and slower/shorter aperiodic fluctuations belong to metallic zigzag and armchair nanotubes, respectively. 



